2005 International Linear Collider Workshop - Stanford, U.S.A. 



Effects of Finite Top Lifetime at the tt Threshold 

Andre H. Hoang and Christoph ReiBer 

Max-Planck-lnstitute for Physics, Munich, Germany 



In this talk electroweak corrections related to the finite top quark lifetime to the total top pair threshold e + e _ 
cross section at NNLL order are discussed. We include the absorptive parts in electroweak matching conditions of 
the NRQCD operators related to the top decay and use the optical theorem. Gauge invariance is maintained. The 
corrections lead to ultraviolet phase space divergences and NLL mixing effects. The corrections can amount to several 
percent and are phenomenologically relevant. 



1. Introduction 

Because in the Standard Model the top quark width T t w 1.5 GeV is much larger than the typical hadronization 
energy Aqcd, the total cross section <r(e + e - — > tt) in the threshold region y/s « 2m t can be computed with 
perturbative methods to very high precision. The rapid rise of the cross section in the threshold region will allow for 
a measurement of the top quark mass (in a threshold mass scheme 1]) with experimental and theoretical uncertainties 
at the level of only 100 MeV. Other parameters such as a 3 , T t or (if the Higgs boson is light) the top Yukawa coupling 
<7tth can also be determined if the normalization and the exact form of the line-shape can be computed with small 
theoretical errors at the percent level Q . This is required because the observable cross section is a convolution of 
the theory prediction with the partly machine-dependent e + e~ luminosity spectrum |2|. 

The theoretical instrument to make first principles predictions for top threshold observables is based on NRQCD, 
an effective theory (EFT) of QCD for heavy quark pairs with small relative velocity v ~ a s <C 1. Within "velocity" 
NRQCD (vNRQCD) which we will use here, it is possible to sum Coulomb singularities and (using renormal- 

ization group methods) logarithms of v to all orders of QCD perturbation theory while a systematic and coherent 
v power counting is maintained. The summation of logarithms avoids large normalization uncertainties 0, E| that 
were obtained in earlier fixed-order predictions 1]. At NNLL order (i.e. including corrections of order v 2 and the 
summation of terms a 2 (a s \nv) n ) all QCD ingredients for the threshold cross section are presently known except 
for still missing subleading mixing effects in the running of the heavy quark pair production current. The current 
normalization QCD uncertainties of the total cross section are estimated to be around 6% [3- 

In this talk we are interested in electroweak effects. At leading order, the three basic electroweak effects are the 
tt production process itself, the finite top lifetime and the luminosity spectrum mentioned above. The latter is 
accounted for in the experimental simulations [2|, the two former effects are described by NRQCD. Here we want to 
discuss finite top lifetime effects at the subleading level including QCD interference effects. Previous partial analyses 
have indicated that the corrections can reach the level of a few percent 0,0. In particular, we investigate the role of 
absorptive parts related to the top quark decay in electroweak loop corrections to the NRQCD matching conditions 
that contribute to the NNLL total cross section. Interestingly, these corrections lead to UV phase space divergences 
that would not exist for stable top quarks and that lead to a NLL renormalization of (e + e~)(e + e~) operators that 
also contributes to the total cross section. For details of our analysis we refer the interested reader to Ref. 0. 



2. Power Counting and Matching Conditions 

Let us first recall the power counting to classify the order at which electroweak effects can contribute by considering 
the matching conditions for the vNRQCD bilinear quark field operators. Electroweak corrections are obtained by 
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matching 2-point functions in the effective theory to those in QCD and the electroweak theory. The NNLL result 
including also all QCD contributions has the form |9| 

£bilinear(*) = j*£>° - (?? + ^ + ^T t (l - - 5m t ^ p {x) + ty p - X p) , (1) 

where the fields i(j p and \v destroy top and antitop quarks with momentum p and positive energies, = (D°, — D) = 
+ igA^ is the ultrasoft gauge covariant derivative and T t is the top quark width defined at the (electroweak) top 
quark pole. The ^-counting is D° ~ m t v 2 ~ I\ ~ m^ 2 since the top width is of order of the typical top kinetic 



energy, which defines the ultrasoft scale mtv 2 . This leads to the counting v ~ cx s ~ g ~ g' for the SU(2) and U(l) 
gauge couplings g and g' , where we can treat the weak mixing being of order one. The term 5m t is a residual mass 
term of order v 2 that arises within threshold mass schemes [lj . Its electroweak contributions are straightforward to 
compute and will not be discussed here further. The term iT t (l — p 2 /2m 2 )) describes the finite lifetime at LL order 
and the NNLL time dilatation effects. It produces the known replacement rule E — > E + iT t to account for finite 
lifetime effects at LL order. Although it renders the effective Lagrangian formally non-hermitian (since hermitian 
conjugation does not change its sign), the EFT inherits unitarity from the underlying theory, so we can later use the 
optical theorem to compute the total rates. It is crucial to understand that the effects of the top decay represent 
hard physics that can be integrated out (i.e. treated by the matching condition iTt) and that, therefore, one can 
determine the corresponding EFT matching conditions for on-shell top quark amplitudes. This is analogous to the 
treatment of photons in an absorptive medium in the optical theory. So to account for all NNLL order electroweak 
effects to Eq. JQ one still has to include the one- loop electroweak and the 0(a s ) and 0(a 2 ) QCD corrections to the 
on-shell top width T t which will also not be discussed here any further. 

Concerning the tt interaction at LL order we only need to consider the Coulomb potential, £ pot = 
J2 PP ' ^-^r^ ^p'^pX-p'X-pi where v is the vNRQCD renormalization scaling parameter. Since effects from 
the top decay are hard (~ m^), we can neglect the momentum exchange p — p' in the electroweak loops during 
the matching computation. From this it is easy to see that all the g 2 (vertex and wave function) corrections to the 
Coulomb potential cancel due to SU(3) gauge invariance. Since all electroweak corrections to the l/m t suppressed 
potentials are beyond NNLL order simply by counting powers of v and g, we also don't have to consider any elec- 
troweak corrections to these potentials. The same argument applies to electroweak corrections to the interactions of 
top quarks with the soft gluons (k ~ m t v) or potential four-quark operators caused by electroweak box diagrams; the 
former can only contribute to 0(a s ) corrections to the potentials, and the latter are already 1/m 2 suppressed without 
even accounting for the powers of g. Moreover, the Coulomb interaction does not generate any UV divergences, so 
electroweak NNLL contributions due to mixing do also not exist. 



The dominant operators used to describe tt spin-triplet production have the form 

®A,p= [e 7,756]^, 



Ov, P = [e 7j e]0' : AiP = [e 7i75 e] OJ, , (2) 



where J p l 



ip P o-j(ia2)X-p - They give the contribution AC = J2 P (CvOv, P + CaOa, p ) + h.c. to the effective 
theory Lagrangian where the hermitian conjugation (which gives the corresponding annihilation operators) is referring 
to the operators, but does not affect their Wilson coefficients. Since we neglect QED radiative corrections, the e ± 
fields act like classic fields, but we need them due to electroweak gauge invariance. The leading order matching 
conditions to Cy/A are obtained from the full theory Born diagrams with photon and Z exchange and are of order 
g 2 . We will see in Sec. |3| that Cy/A will receive imaginary matching conditions from one- loop electroweak corrections 
that contribute at NNLL order and play a role similar to the imaginary terms in Eq. Q}. 

It has been shown in Refs. that up to NLL order for a tot there are no QCD interference effects from (ultrasoft) 
gluon radiation off the top/antitop quark or its decay products. (For this consideration only ultrasoft gluons are 
relevant because they cannot kick a top quark off-shell.) The proof was conducted by explicit computation of 
diagrams. Using the analysis of matching conditions we can show that this statement even holds at NNLL order. For 
the time-like ultrasoft A gluons, QCD gauge invariance ensures that the dominant electroweak matching corrections 
to the A interaction vertex shown in Eq. vanish because we can again approximate the ultrasoft gluon momentum 
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Figure 1: (a) Full theory diagrams in Feynman gauge needed to determine the electroweak absorptive parts in the Wilson 
coefficients Cy/A related to the physical bW + and bW~ intermediate states. Only the bW + cut is drawn explicitly, (b) Full 
theory diagrams describing the process e + e~ — > bW + bW~ with one or two intermediate top or antitop quark propagators. 
The circle in the first diagram represents the QCD form factors for the tt vector/axial- vector currents. 

in the electroweak loop as zero. Moreover, the exchange of time-like ultrasoft A gluons does itself not contribute 
at LL because they can be removed from the LL particle-antiparticle sector in Eq. (Q by a redefinition of the top 
and antitop fields related to static Wilson lines 13]. An anomalous interaction in analogy to the g — 2 is not covered 
by this argument but suppressed by a factor 1/raf. Accounting also for powers of g one finds that interference from 
A gluons does not contribute at NNLL order. For the space-like ultrasoft A gluons the conclusion is the same 
because they couple to the quarks with the p.A/m t coupling. The g 2 -suppression from an additional electroweak 
loop correction to the interaction vertex then leads to a contribution beyond the NNLL order. 



3. Absorptive Matching Conditions 



To determine the absorptive parts related to the top decay of the matching conditions of the operators Oy, P and 
Oa, p that contribute to the total cross section at NNLL order, we have to consider the bW + and bW~ cuts of the 
full theory diagrams shown in Fig.^Ji. To obtain the contributions at this order the external (on-shell) top quarks 
can be taken to be at rest. The full theory amplitude has the form 



A 



v e+ (k')^(iC] 



bW,abs 
V 



CbW,abs \ 
A 7s) U e 



(k) u t {p)-f^vt(p) 



(3) 



where k + k' = 2p = (2m t , 0). The amplitude of the charge conjugated process describing top pair annihilation reads 



A 



u e -(k)r(< w ^ 



+ iC A ' a s 75) v e + (k 1 ) Vt(p) 7 M u t (p) 



(4) 



We used the cutting equations to obtain expressions for zCj^' abs and checked electroweak gauge invariance by 
carrying out the computation in unitary and Feynman gauge. Analytic formulae are provided in Ref. 10]. The 
results are consistent with results obtained earlier in Ref. |8|. It is a consequence of the unitarity of the underlying 
theory that the sign of the imaginary part of the amplitude does not change in the charge conjugated amplitude. It 
is straightforward to match the amplitudes for the operators Ov/a, p an d Oy/A P ^° ^ ne ^ un theory results in Eqs. (|3J) 
and (El). The resulting matching conditions for the operators Ov/a n and 0\ rl A read 



\. The resulting matching conditions for the operators Oy/A,p an d ^y/A p 

C v/A {y=l) 



z^born | • ^bW,abs 

"I" 1 ^v/a • 



(5) 



where we have also indicated the Born level contributions. In a full treatment of electroweak effects the coefficients 
Cy/A also include the real parts of the full set of electroweak one- loop diagrams indicated in Fig.[IJi. A comprehensive 
examination of these contributions will be provided in subsequent analyses. 



0402 



4. Time-Ordered Product and Renormalization 



Through the optical theorem the NNLL order corrections to the total cross section that come from the absorptive 
one-loop electroweak matching conditions for the operators Oy/A,p an d from the time dilatation corrections can be 
computed from the imaginary part of the (e + e~)(e + e~) forward scattering amplitude, 



Otot 



1 



Im 



(C^) + C») L lk A[ k 



where L lk = \ (k + k') 2 (S lk — e l e k ) is the spin-averaged lepton tensor and 



A[ k =iJ2 \ d ^ x e ~ i( *' X ( | T ° l vl (°) °pM ( x ) I ) = 2 7V C G°(a, v, m u v) 



(6) 



(7) 



,7 t 



is the time-ordered product of the tt production and annihilation operators 3 p l and J p l _6]. Here we used k + k' = 
(-v/s.O), e = k/|k| and q = {\fs — 2mt,0), \fs being the cm. energy. The result reads 

A<r t r i = 2Ar c Im|2i[Cr n ^ W ' abs + ^ orn ^ W ' abS ]G (a,«,mi,i,)+ [ (C£ ora ) 2 + (C^") 2 ] SG$(a, v, m t , v) J , (8) 

where a = — Vc^(^)/47r = Cpc^^t^)- The term G° is the zero-distance S-wave Green function of the non-relativistic 
Schrodinger equation which is obtained from the LL terms in Eqs. © and the Coulomb potential. In dimensional 
regularization it has the form [3 



G°(a,v,m t ,v) 



m 
4tt 



In 



-+ln2 + 7£ + V> 1-- 



mf a 1 
47 



47T 



(9) 



where v = ((v^ — 2(mt + Jmt) + irtj/mt) 1 ^ 2 . The term JGp represents the corrections originating from the time 
dilatation correction in Eq. JU and reads £Gp(a, mt, v) — — + f ^7 + a ^;]G (a, m^, Note that the 

Wilson coefficients do not have a LL anomalous dimension, so only the matching conditions at v — 1 appear 

in Eq.©. 

One can check that the terms proportional to Cy^' abs in Eq. © are in agreement with the full theory matrix 
elements from the interference between the double-resonant amplitudes for the process e + e~ — > tt — > 
(first diagram in Fig.^J)) and the single-resonant amplitudes describing the processes e + e~ — ► t + — > bW + bW~ 
and e + e~ — > bW + bW~ in the threshold limit for m t 00 (subsequent diagrams in FigJIJ)). To find 

literal agreement between full and effective theory matrix elements one has to replace the ie terms in the resonant full 
theory top propagators by the Breit-Wigner term im t T t /2. As discussed above, there are no further QCD corrections 
in the non-relativistic limit due to the cancellation of the QCD interference effects caused by gluons with ultrasoft 
momenta. 

The corrections given in Eq. © exhibit UV 1/e-divergences that have interesting features. Physically they arise 
from a logarithmic high energy behavior of the top-antitop effective theory phase space integration for matrix elements 
containing a single insertion of the Coulomb potential. Technically they enter the imaginary part of the forward 
scattering amplitude because the imaginary parts of the matching conditions of Eqs. ((5J) lead to a dependence on the 
real part of G° (see Eq. ©). In the full theory this logarithmic behavior is regularized by the top quark mass. It 
is important that the divergences only exist because the top quark is not treated as a stable particle. In particular, 
the UV divergences from the time dilatation corrections arise from the Breit-Wigner-type high energy behavior of 
the effective theory top propagator derived from Eq. ([T]). which differs from the one for a stable particle. Likewise, 
the interference effects described by the absorptive electroweak matching conditions for the operators Oy, P and Oa, p 
would not have to be accounted for if the top quarks were stable particles. UV divergences of the same kind for 
the NNLL total cross section have been observed and noted before fl B Q Q, but no resolution of the issue was 
provided. From the point of view of having an EFT with non-hermitian contributions it is obvious that these UV 
divergences must be handled in the canonical way using renormalization. Since the divergences are directly related 
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to operators with non-hermitian Wilson coefficients, the renormalization procedure will naturally involve operators 
with non-hermitian Wilson coefficients. 

The operators that are renormalized by the UV divergences displayed in Eq. (0) are the two (e + e~)(e + e~) operators 

O v = ~ [e^e] [e7 M e] , & A = - [e 7 M 75 e ] [e 7/i 7 5 e] , (10) 

which give the additional contribution AC = CyOy + CaOa to the effective theory Lagrangian, Cy/A being the 
Wilson coefficients. Because in this work we neglect QED effects, the electron and positron act as classic fields and 
therefore Cy and Ca run only through mixing due to UV divergences such as in Eq. @ . Since only the imaginary 
parts of the coefficients Cy/A are relevant for the discussion, we neglect the real contributions in the following. Using 
the standard MS subtraction procedure to determine the (non-hermitian) counterterms of the renormalized Oy/A 
operators and standard methods to compute and to solve the anomalous dimensions, one obtains the following form 
of the Wilson coefficients Cy/A fc> r scales below mt (v < 1), 



Cy/ A (y) = c v/A (i) + i c * 



I V^v/A ) y^y/A ) ) ~~ "I" ^y/A 

\ / lilt 



-,bW,abs 



m, ' - y l A 



ln(z) 



^ C F r t / r borm2i 2/^ , A + 2C F ) ?t /^born^ / n 1 -.n 

Po m t ' m t 1 ) 

where z = a s (m t iy)/a s (m t ), p(z) = ^ — ^ In 2 2 + ln21n(z) — Li 2 (|) and the CV/a(1) are the hard matching 
conditions. The operators OyM l ea d to an additional contribution to the total cross section of the form 



Acr t F ot = Im 



Cy + 



(12) 



For details of the computations see Ref. Parametrically Acr^f is of order g 6 and thus constitutes a NLL 

contribution as one can also see from the fact that the corresponding UV divergences were generated in NNLL 

r.2 

7 tot 



order effective theory matrix elements. The correction Aa[ ' t 2 is energy-independent, but it is scale-dependent and 



compensates the logarithmic scale-dependence in the NNLL contribution Aa^ 1 . The matching conditions Cv/a( 
1) are presently unknown and will be analysed in subsequent work. For now we set them to zero in the numerical 
analysis presented below. We note that it was shown in 9] that the difference between the full theory phase space 
(which is cut off by the large, but finite m t ) and the effective theory phase space (which is infinite in the computation 
of the forward scattering amplitude) contributes to Cy/A^y — 1) and also represents a NLL effect. 



v 



5. Numerical Analysis 

In Fig.[2]we have plotted the sum of Ad^ 1 and Acrjj t 2 in picobarn in the IS mass scheme |£j for Mis = 175 GeV, 
the fine structure constant a — 1/125.7, s 2 w — 0.23120, Vtb = 1 and Mw = 80.425 GeV with the renormalization 
scaling parameter v = 0.1 (solid curves), 0.2 (dashed curves) and 0.3 (dotted curves). For the QCD coupling we 
used a s (Mz) = 0.118 as an input and employed 4-loop renormalization group running. Note that in the IS scheme 
Suit = Mis(Vc^(^)/47r) 2 /8. For the top quark width we adopted the value T t = 1.43 GeV. Note that in a complete 
analysis of electroweak effects the top quark width depends on the input parameters given above and is not an 
independent parameter. For the purpose of the numerical analysis in this work, however, our treatment is sufficient. 
We find that the sum of the corrections is negative and shows a moderate v dependence. We find that the corrections 
are around —10% for energies below the peak, between —2% and —4% close to the peak and about —2% above the 
peak. Interestingly, they partly compensate the sizeable positive QCD corrections found in 0,^3- The peculiar 
energy dependence of the corrections, caused by the dependence on the real part of the Green function G°, also leads 
to a slight displacement of the peak position. Relative to the peak position of the LL cross section one obtains a 
shift of (30,35,47) MeV for v = (0.1,0.2,0.3). This shift is comparable to the expected experimental uncertainties 
of the top mass measurements from the threshold scan Q. 
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Figure 2: The sum Aa^ 1 + Aa t r ' t 2 in pb for Mis = 175 GeV, a = 1/125.7, s 2 w = 0.23120, V tb = 1, M w = 80.425 GeV, 
Ft = 1.43 GeV and v — 0.1 (solid curves), 0.2 (dashed curves) and 0.3 (dotted curves) in the energy range 346 GeV < y/s < 
354 GeV. 
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